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Abstract 

Nested canalizing Boolean functions play an important role in biological motivated 
regulative networks but also in signal processing, such as in describing stack filters. 
It has been conjectured that this class of functions has a stabilizing effect on the 
network dynamics. It is well known that the average sensitivity plays a central 
role for the stability of (Random) Boolean networks. Here we prove a tight upper 
bound on the average sensitivity for nested canalizing functions in dependence of 
the number of relevant input variables. We further show that it is smaller than | as 
conjectured in literature. This shows that a large number of functions appearing in 
biological networks belong to a class that has a very low average sensitivity, which 
is even close to a tight lower bound. 

1. Introduction 

Boolean networks play an important role in modeling and understanding signal trans- 
duction and regulatory networks. These networks have been widely studied with focus 
on many facets, e.g. [TJ EJ [3]. One line of research focuses on the dynamical stability of 
randomly created networks. For example, random Boolean networks tend to be unstable 
if the functions are chosen from all possible Boolean functions and the average number 
of variables (average in-degree) is larger than two This can be attributed to the 
fact that expected average sensitivity of the random functions, which is an appropriate 
measure for the stability of random Boolean networks 0(6], is too large. 
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If only functions from certain classes are chosen, stable behaviour can be achieved for 
higher in-degrees. For instance, canalizing and nested canalizing functions, introduced 
in [7J [8], have been conjectured [9] to have a stabilizing effect on network dynamics. In 
|10j it has been shown that stable behaviour is possible even for networks with in-degree. 
Interestingly, studies of regulatory network models have shown that a large number of 
their functions are canalizing [Til 021 031 ttH US US]- Canalizing functions are also 
important for the construction of stack filters used in signal processing |17j . 

A function is canalizing in a variable, if its output is constant when this variable is 
set to its canalizing value. Nested canalizing functions are canalizing functions, whose 
restriction to the non-canalizing value is again a canalizing function and so on (a precise 
definition is given later). In this paper we analyze nested canalizing functions and 
we focus in particular on their average sensitivity. The notion of sensitivity was first 
introduced by Cook et al. [18J. Later it was applied to Boolean functions [19] and can be 
viewed as a metric for the influence of a random permutation of the input variables on the 
output of the function. Since then the average sensitivity has been studied widely. In [20] 
the average sensitivity in the context of monotone Boolean functions was investigated 
and an upper bound for locally monotone functions was presented in [15]. Here we give a 
tight upper-bound on the average sensitivity of nested canalizing functions with different 
numbers of variables. This also shows that the average sensitivity is always smaller than 
| as conjectured in J2JJ. We further give a recursive expression of the average sensitivity 
and the bias of this class of functions. Finally we will discuss and compare our new 
bounds and some old bounds. 

Our main tool is the Fourier analysis [221 123] of Boolean functions, which is introduced 
in Section [2j There we also address further concepts needed, such as restrictions of 
Boolean functions. In Section [3] spectral properties of canalizing and nested canalizing 
functions are broached. Additionally we discuss functions, in which all variables are most 
dominant, as they turn out later to minimize the average sensitivity. In Section U] the new 
bounds on the average sensitivity are presented based on a recursive expression of the 
average sensitivity of nested canalizing functions. We conclude then with a discussion 
of the obtained results and some final remarks. 

2. Notation, Basic Definitions and Fourier Analysis of Boolean 
Functions 

A Boolean function (BF) / € T n = {/ : VL n — > £2} with Q = {— 1,+1} maps n-ary 
input tuples to a binary output. In general not all input variables have an impact on 
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the output, i.e. are relevant. 

Definition 1. ]21^ A variable i is relevant to a BF f , if there exists an x 6 £l n such 
that 



where x © a is the vector obtained from x by flipping its i-th entry. 
Further rel(f) is the set containing all relevant variables of f. 

2.1. Fourier Analysis of Boolean Functions 

In this section we will recall some basic concepts of Fourier analysis and some results 
concerning restrictions of BF as shown in [Mj. Let us consider x = (x\, X2, ■ ■ ■ , x n ) 
as an instance of a uniform distributed random vector X = (Xi,X2, ■ ■ ■ ,X n ), i.e., its 
probability density functions can be written as 



We first define the basis of the Fourier transform as the set of monomials xu( x ) given 
by all U C [n] = {1, 2, . . . , n}, where 



For U = we set X0( x ) = 1- 

It is well known that any BF / can be expressed by the following sum, called Fourier- 
expansion [22l [23] , 



/(x)^/(x© ei ), 



Pr[X = xl = — . 



(1) 



/(x)= W)-Xu(x.) 



UC[n] 



where f(U) are the Fourier coefficients. The Fourier coefficients can be recovered by 




(2) 



X 



Let A C U and A = U\A, then 



Xu(x) = Xa(x) • Xa( x ) 



which directly follows from the definition of xu (Eq. ©)• 
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2.2. Restrictions of Boolean Functions 

If we restrict /, i.e. if we set the i-th input variable of / to some constant a £ {—1, +1}, 
we express the obtained new function by /(*' a ) £ T n . Every BF can be decomposed 
in two unique restricted functions for each relevant variable, as stated in the following 
proposition: 

Proposition 1. For any f £ T n and each i £ [n] there exist unique functions f^' + \ f^ 1 '^ £ 
T n , with i ^ rel(f( l ' + ^) and i ^ rel(f^ l ~^), such that 

f = g (i,+) f (i,+) +g (i,-) f (i,-) j 
where the functions g^' + \ g^'^ £ T n are given by 




(*>+)(x) = 1 i and /- } (x) 



1 if Xi 
else 



The Fourier coefficients of these restricted functions can be derived as stated in the 
following proposition, the proof and an extension for product distributed input variables 
can be found in 



Proposition 2. \2$ Let f be a BF in n uniformly distributed variables. Consider the 
restricted function obtained by setting Xi = a\, then 

p«)(U) = j(U) + a i -f(UU{i}) 

where U C [n] \ {i}. 

The reverse relation, i.e. the composition of a BF by two restricted functions, is 
described in terms of Fourier coefficients by the following proposition. The proof can 
again be found in 



Proposition 3. fE$ The Fourier coefficients of any BF f with uniform distributed input 
variables can be composed of the coefficients of its two restricted functions and 
/(*.+!) 

/(tO = \ (f (t ' +) (u \ {%}) + (-i)i^»i/(*,-) (c/ \ W) ) , 



or 

' P' + \U\{i}) + p>-\U\{i}) ifi£U 
fd,+)(y)- p>-)(U) ifi$U 



2/(17) 
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The zero coefficient /(0) plays an import role in the analysis of BFs. In the uniform 
case it corresponds to the bias of the function /, where bias is defined as the probability 
that f{x) is —1. Next, we show that it can be easily composed by the zero coefficients of 
the restricted functions. We will need this later to show the recursive behavior of nested 
canalizing functions. 

Corollary 1. The zero coefficient of any Boolean function f can be written as: 

fm = \f {i ' +) m + lf {i '- ] m, (3) 

where i G [n] is the index of some variable. 

Proof. Follows directly from Proposition [3l □ 

If we restrict a function to more than one variable, namely to a set of variables K, we 
denote the restricted function with f^ K ' a \ where a is a vector containing the values to 
which the functions is restricted. The Fourier coefficients are then given by the following 
proposition, which again can be found in [24]. 

Proposition 4. [2^] Let f be a Boolean function and f(U) its Fourier coefficients. 
Furthermore, let K be a set containing the indices i of the input variables Xj, which are 
fixed to certain values aj. The Fourier coefficients of the restricted function are then 
given as: 

fi^)(U)= (* s (a)./(tfU5)), 

SCK 

where a is a vector containing all ai,i E K. 

3. Nested Canalizing Functions 

3.1. General 

To define nested canalizing functions (NCF) we first need to look at canalizing functions: 

Definition 2. A BF f is called < i : a : b > canalizing if there exists a canalizing 
variable Xi and a Boolean value a £ { — 1,+1} such that the function 

/(x| Ii=fl ) = /M( x) = 6i (4) 
for all xi, ...Xi-i, Xi + \....x n , where b € {—1, +1} is a constant. 
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Hence, / is canalizing in variable % if and only if fractionation according to Proposition 
[1] results in either f^' + ^ or f^ 1 '^ being a constant function. 

As shown in [23J the Fourier coefficients then fulfill the following condition: 

/(0) + o i ./({i})=6 i . (5) 

A NCF can be described recursively as a canalizing function, whose restriction is again 
nested canalizing or more formally: 

Definition 3. For k = 1 and k = any BF with k < n relevant variables is a NCF. 
For k > 1 a BF is a NCF iff there exists at least one variable i and two constants 
cti,Pi e {+1,-1} such that = ft and /(»>-"*) is a NCF with k - 1 relevant 

variables. 

Further let aWi); • • • > x ir(k) De the variable order for which a NCF fulfills the properties 
from this definition, then we call, according to [21] . such a function {tt : a : f3} nested 
canalizing. 

This condition can be easily expressed using the Fourier representations: |24j 
/ is {tt : a : f3} nested canalizing, if for all j G {1, . . . , k} 

E K n{i} ' • Xs\ W (a) ■ /(<§)) = fa 
SQj] 

where a is a vector containing all negated aj, i.e. on = —cti and S is a set which is 
retrieved by applying the permutation tt to the elements of S. 

Further we may want to add an example: Let / be {tt : a : f3} NCF with k = 2 
relevant variables and tt = id, that is S = tt(S) = S, then we can write: 

ft = /(0) + «i/({i}) 

ft = /(0) - aj({l}) + a 2 /({2}) - aia 2 /(l, 2). 
3.2. Properties of Nested Canalizing Functions 

In this section we state some properties of NCF. First we address most dominant vari- 
ables, which are defined as: 

Definition 4. According to 121)/ we call a variable i a most dominant variable of f if 
there exists at least one variable order tt = («,...), for which f is {tt : a : ft} canalizing. 

The set of most dominant variables has an impact on a number of Fourier coefficients, 
which is summarized in the following proposition. 
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Proposition 5. Let K be the set of most dominant variables of a {it : a : /?} NCF f 

with uniform distributed input variables. Then the absolute values of the corresponding 
Fourier coefficients are all equal the same constant c > 0, i.e., 



f(U) 



VUQK,U^ {0}, 



or, more general, 



a 



3 ■ Xu\{ j} (a) ■ f(U) =c VC7 C K,U + {0} and Vj G K. 



(6) 



Further the absolute value of the zero coefficient, /(0), is 1 — c and the sign is given by 
b, i.e., 

b = sgn (/(0) N 
and 

f3i = b VieK. 
The proof can be found in Appendix 1X1 

For the special case, in which all variables are most dominant, we derive the following 
two corollaries: 

Corollary 2. Let f be a {ir : a : /3} NCF with n uniform distributed and k relevant 
input variables. All variables are most canalizing if and only if the Fourier coefficients 
fulfill the following conditions, 

ctj ■ j Yl «i J f( S ) = c V5 C [n], 5 ^ {0} and Vj G 5 (7) 

and hence 



f(S) =c VSCK,S^ {0}, 



7(0) 



(8) 



1 - c 



with 



c = 2- {k - l \ (9) 

Proof. Eq. ([7]) and ([Sj) follow directly from Proposition [5l while Eq. (|9|) follows from 
Parsevals theorem. 
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□ 



Corollary 3. Let f be a {ir : a : /?} NCF with k > 1 uniform distributed and relevant 
input variables. All variables are most canalizing and fa = b, \/i € {1, . . . , k}. All such 
NCFs are completely described by a and b. and hence there are 2^ k+1 ^ such functions. 

Proof. The proof follows directly from the previous Corollary. □ 

As mentioned before, the zero coefficient plays an important role. Interestingly, we 
can describe the zero coefficients for NCFs in a recursive manner: 

Proposition 6. The zero coefficient of a {ir : a : /?} NFC f can be recursively written 
as : 

; (0) = i />*.«*) (0) + lfa. 

Proof. Follows directly from Corollary [TJ □ 

Further, the zero coefficient is upper-bounded as shown in the following proposition: 

Proposition 7. The absolute value of zero coefficient of a {tt : a : (3} NCF f with 
uniform distributed input variables can be bounded as: 

^(^(-l) fe + l) < 1/(0)1 <1-2^T> 
where k = rel(f) is the number of relevant variables. 

Proof. First, we prove the right hand side: Using the triangle inequality we get from 
Proposition 

|/(0)| <i|/>««)(0)| + i. 

Obviously the zero coefficient of a function with only one relevant variable % is zero. The 
proposition now follows using induction. The left hand side can be easily shown using 
the inverse triangle inequality and induction. □ 

As seen in Corollary [2J a NCF, whose variables are most dominant, fulfills the upper- 
bound with equality. Further, it can be easily seen that NCFs with alternating fa, i.e., 
with (3 = (-1, +1, -1, +1, . . .) or f3 = (+1, -1, +1, -1, . . .) fulfill the lower-bound with 
equality. 
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4. Average Sensitivity 



4.1. Definition 

The average sensitivity (as) [19] is a measure to quantify the influence of random per- 
turbations of the inputs of Boolean functions. It is defined as the sum of the influences 
of the inputs of the function, which is defined as the probability of a change of the 
function's output if input i is flipped: 

Definition 5. f}25\ \26 $ ) Define the influence of variable i on the function f as 

I i (f) = Pr[f(X)^f(X®e i )}. 
The influence can be related to the Fourier spectra as follows [27]: 



SQnhieS 



The average sensitivity is defined as the sum of the influences of all input variables of 
Definition 6. (125\, \19jj ) The average sensitivity of f to all input variables is defined as 

35(f) = £ h{f). 

ie[n] 

Consequently the average sensitivity can also be expressed in terms of the Fourier 
coefficients [25] : 

as(/)= £ f(S) 2 \S\. (10) 

4.2. Restricted Functions 

To investigate the behavior of the average sensitivity of restricted functions we first need 
to define the function £ : T n X T n — > R by 



ZV,g) = \ fi- E hu)S(u)] ■ 

\ UC[n] J 



We can then state the following theorem, which shows the relation between the average 
sensitivity of a BF and the average sensitivity of its two restricted functions. 
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Theorem 1. Let f^ l ' + \f^ 1 ' be the restrictions of f to some relevant variable i of f. 
Then 

as(f) =\as{f^) + \as{f^) +£(/(*>+),/(<>-)) 

The proof can be found in Appendix [Bj For NCFs we obtain then: 

Corollary 4. The average sensitivity of a {it : a : /3} NCF can be recursively described 
as: 

as(/) = \ (as(/^«)) + 1 - />* a <>(0)/3i) . (12) 

In [21] an upper-bound on the average sensitivity of NCF has been conjectured. In 
the following theorem, we proof this conjecture to be correct. 

Theorem 2. The average sensitivity of a NCF with k = rel(f) relevant and uniform 
distributed variables is bounded by 

~^J~i < «(/) < | " 2"* " | • 2- fc • (-l) fc . (13) 

The bounds are tight. 

The proof can be found in Appendix O The tightness is shown in Corollaries [6] and [71 
We can generalize the upper bound of Theorem [2] as: 

Corollary 5. The average sensitivity of a NCF with k = rel(f) relevant and uniform 
distributed variables is upper-bounded by 

-</) < f . 

Corollary 6. Let f be a NCF, whose variables are all most dominant, then f fulfills 
the left-hand-side of Theorem with equality. 

Proof. We start from Corollary [H using the fact that in this case |/(0)| = 1 — t^t an d 
all ^ = sgn(f(%)), we get: 

as(/) = \ (as(/(^)) + 1 - (1 - -^)) (14) 

= \ (as(/^) + -^)) • (15) 
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Since as(/) has k relevant variables, while as(/fa' Qi )) depends only on k — 1 relevant 
variables, we can also express Eq. (fHj) as: 

as(fc) = 1 (as(fc - 1) + JL)) . 

The proof can now be concluded by solving this recursion using induction. □ 

Corollary 7. Let f be a NCF with alternating ft, i.e. f3 = (— 1, +1, — 1, +1, . . .) or 
(3 = (+1,-1, +1,-1,...). Then f fulfills the right-hand-side of Theorem^ with equality. 

Proof. According to the proof of the previous corollary we start from Corollary 0] and 
use |/(0)| = | (-^hri— l) fe + l)j the poof is concluded by solving the recursion. □ 

The last two corollaries showed that the maximal average sensitivity is achieved if the 
bias, i.e. the zero coefficient, is minimized, and vice versa. Next, we derive and common 
bound as given in the following proposition: 

Proposition 8. Let f be any NCF and /(0) its zero coefficient, then 

as(/) + |/(0)| < |. 
Proof. Combining Corollaries H] and we get: 

aS( /(^))_ ft ./(^)(0)<| 

and since ft G { — 1,+1}: 

as (/(^)) + l/fa- 1 *) (0)| < -. 

3 

Substituting /fa> a *) by / concludes the proof. □ 

5. Discussion 

In Figure [T] we summarized the most important bounds from the previous section. We 
plotted the average sensitivity versus the bias. Additionally we included a general lower 
bound on the average sensitivity as it can be found in |15j . One can see that this bounds 
intersects with our lower bound (which we plotted for k = 5), though we stated that our 
bound is tight. However, this is not a contradiction, since the lower bound of Theorem 
[2] is only achieved for highly biased functions, which are located outside the intersection. 
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CorollaryE] ----Theorem El k=5 (LHS) 

lower bound [15 - - Proposition [8] 



1.5 



^ 1 
0.5 


-1 -0.5 0.5 1 

Z(0) 

Figure 1: Important Bounds: The dotted-area gives the possible values for the average 
sensitivity of any NCF, the lined area for BF with k = 5 input variables. 

For k = 5 our lower bound forms a triangle with the upper bound as formulated in 
Proposition [8l where the NCF with all variables being most dominant are located in the 
left and right corners, however, for larger k the lower bound will move further to the 
bottom, and with them the most dominant NCFs. 

The upper bound from Corollary [5] also intersects with the bound formulated in Propo- 
sition [HJ Again, this is not a contradiction, since NCF reach this bound only for small 
bias. 

In general the average sensitivity is upper-bounded by k, i.e., as(/) < k. Fur- 
ther, as shown in [15], the upper-bound for unate, i.e., locally monotone, functions 
is as(/) < \J (1 — f(%))k. A function is a member of the class of unate functions, if it is 
monotone in each variable. That means in terms of regulatory networks that an input to 
a function, i.e., a regulator, can be either an activator or an inhibitor towards a certain 
gene. Further, NCFs form a subclass of the unate functions. Thus, the average sensitiv- 
ity of NCFs is remarkably low. Since a low average sensitivity has a positive influence 
on the stability of Boolean networks [2], we can conjecture that networks consisting of 
NCF are more stable. 
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6. Conclusion 



In this paper we investigated Boolean functions, in particular canalizing and nested 
canalizing functions, using Fourier analysis. We gave recursive representations for the 
zero coefficient and the average sensitivity based on the concept of restricted BFs. 

Further, we addressed the average sensitivity for nested canalizing functions and de- 
rived and proofed a tight upper and lower bound. We show that the lower bound is 
achieved by functions, whose input variables are all most dominant and which are maxi- 
mizing the zero coefficient and, hence, the bias. The upper bound is reached by functions, 
whose canalized values are alternating, and which are minimizing the bias. 

We then generalized the upper bound to as(/) < |, which has been conjectured in 
literature, however no proof has been given so far. Finally we derived a common bound 
for bias and average sensitivity and discussed the stabilizing influence of the class of 
nested canalizing functions on the network dynamics. 

It is worth noting that all these results depend on the assumption of uniform dis- 
tributed inputs (see Section [2|). This opens the question if the results can be generalized 
to other distributions. The recursive representations can be easily extended to prod- 
uct distributed input variables. But without further constraints there always exists a 
distribution which maximizes the (accordingly defined) average sensitivity, i.e., for any 
function with k relevant variables the average sensitivity can be k. 
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A. Proof of Proposition [5] 

Proof. The proof for the zero and first order coefficients, i.e. \U\ = 1 and U = 0, 
follows directly from Eq. ([5]). We can hence use Eq. (J6j) as an induction hypothesis for 
coefficients with order smaller than \U\. We show next that as a result this is also valid 
for coefficients with order \U\ + 1. 
Since 

f(UA) {T)= Y j (»(«)-/(TuS)), 

scu 
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and / is canalizing in any variable k, and hence every restriction of / must also be 
canalizing in variable k, i.e. /^(0) + a k • f^ U \{k}) = b, we get: 

b = Yl {xs(a) ■ f(S)) +a k -Y, (xs(«) ■ /(« U 5)) , 
scu scu 

b = E (-^•x5\m(«)-/(s))+/(0) 

SCU,S^$ 

+ a k ■ Yl (xs(a) ■ f({k} U S)) + a k ■ /({*}). 

scu,s^% 

Using Eq. (JSJ) and using the induction hypothesis, we get 

°= E (" c )+ E («*-Xs(«)-/(Wu5)) 
+ (a* • Xc/(«) • f({k} U 17)) . 
We again assume (0) as true for all S C U, i.e. |5| < |J7| and, hence, write: 

= _( 2 M - 1) • c + (2^1 - 2) • c + (a* • X c/(a) • /({&} U 17)) 
c = a k ■ Xuipc) ■ f({k} U U), 

which concludes the proof. □ 

B. Proof of Theorem Q] 

Proof. Starting from the definition of as as given in Eq. (|10p . we can fractionize the 
Fourier coefficients according to Proposition This yields in: 

as(/) = E (\f il ' +) ( S \ ») + 5 (-l) |SnWI /^(S \ «)) 2 

=\ E (((+D |5n{i}l / (i ' +) (^\w)) 2 

+ ((-l) |5n{i}l f l '-\S \ {i})) 2 + 2 (-i)l* n WI /(*.+) (5 \ {*})/<*'-> (S \ {»})) |5| 
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which leads us to: 



as ^)=I E (f (i ' + \s\{i})) 2 \s\ 

SC[n],SyfD 

+\ E {f (i '~\s\{i})Y\s\ 

SC[n],S^9 

+\ E (-if mi f {i ' + \s\{i})p>-\s\{i})\si 

SC[n],S^$ 



and hence to: 



(f)=\ E {f^ +) (s)) 2 \s\ + \ Yl (f {t ' +) (S))\l + \S\) 

SC[n]\{i},S^ SC[n]\{i} 

+\ E (f {i '-\s)) 2 \s\ + \ Yl (f {l '~ ) (S))\l + \S\) 



SC[n]\{i},S^Hl SC[n]\{i} 

+\ £ (-i)°P' + \s)f^-\s)\s\ 

SC[n\\{i},S^9 

4 E f^ + \S)p>-\S) (1 + |5|). 

SC[n]\« 

Since / (i,a) (S) = for all S" : i G S 1 we can write 

E (/ (i - +) (^))Vi+^ E 

SC[n],5^0 5C[n] 

<> .. ' <> 



(/«,+>) =1 

2 



+^ E (/^(^Vi+^E 

5C[n],5^0 SC[n] 



=as(/«.-)) =1 



SC[n],S^0 SC[n],S^0 

~ E f^wP-Hs)- 

SC[n] 



Finally we get 



as(/) 4as(/(*' + )) + ^as(/^-)) + 1 - \ £ P> + \S)f^~\S), 



2 yJ ' 2 VJ y 2 2 

5C[n] 
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which concludes the proof. □ 

C. Proof of Theorem [2] 

Proof. We first prove the right hand side. Let us recall Corollary UJ 

as(/) = \ (as(/^«)) + 1 - f^Xm) ■ 
If we apply it again and use Proposition [6j we get: 

as(/) =\ Q ( as (/(^) (7r * +1 ' ai+l) ) + 1 - /(W^^ (0)A+1 ) 

= ias(/^^) K+1 ' Qs+l) ) + - - J/>^) (,ri+1 ' ai+l) (0)ft +1 

= l as(/ (^ Ql )(- + — )) _ l/tW-m.-* i) (0)(A + 
1 3 

- /mft + 4 

Since ft, ft+i G {-l,+l}and|/>-^) K+1 ' Ql+l) (0)| < 1, we can upper-bound _ l/>i.«i) ( " i + 1 ' a<+l) (0)(ft + 

ft+i) - as: 

-^/ ( — )( " +1 '" +l) (0)(ft + ft+i) - < J 

and finally upper-bound as(/) as 

as(/)<ias(/^«) (7r<+1 ^ +l) ) + l, (16) 

where y ? ( 7r »> Q! ») 7r * +1,Wt+1 has k — 2 relevant variables. We will now show the theorem by 
induction. For k = 1 Eq. ()13j) simplifies to 

as(/) < 1, 

which is obviously true. For k = 2 Eq. (fT3|) results in 

as(/) < 1, 
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which is also true and can be verified by inspecting all possible functions. 

Using Eq. (|13j) as the induction hypothesis, we will now show that Eq. 
for k if it is true for k — 2. Using Eq. [TBI we get: 



CED is true 



as(/)<^s(/^^ lT ' +1 ' a ' + Vl 
since y( 7r i> a i)^ 7ri+1 ' a ' +1 ' 1 h as k — 2 relevant variables, we can use our hypothesis and write 

as(/) < \ (| - 2-( fc - 2 > - \ ■ 2-( fe - 2 ) • (-1)*' 2 ) + 1 

= - -2~ k --2- k (-l) k , 
3 3 V ' ' 

which concludes the induction. 

The left hand side is commonly known in literature and can be proven along the lines 
like the right hand side using the following inequality, which follows from Corollary [J] 
and Proposition [7J 

as(/) > \ (as(/(— )) + ^) • 

□ 

References 

[1] M. Aldana, Boolean Dynamics of Networks with scale- free Topol- 
ogy, Physica D: Nonlinear Phenomena 185 (1) (2003) 45-66. 
|doi : 10 . 1016/S0167-2789 (03) 00174-X[ 

[2] I. Shmulevich, S. A. Kauffman, Activities and sensitivities in boolean 
network models, Physical Review Letters 93 (4) (2004) 048701. 
|doi : 10 . 1 103/PhysRevLett . 93 . 4870 1 . 

[3] S. Bornholdt, Boolean network models of cellular regulation: prospects and limita- 
tions, Journal of the Royal Society, Interface / the Royal Society 5 Suppl 1 (2008) 
S85-94, PMID: 18508746. |doi : 10 . 1098/rsif . 2008 . 0132 . f ocusl 

[4] S. A. Kauffman, Metabolic stability and epigenesis in randomly constructed nets, 
Journal of Theoretical Biology 22 (1969) 437-467. 

[5] J. F. Lynch, Dynamics of random boolean networks, in: R. C. K. Mahdavi, 
J. Boucher (Eds.), Current Developments in Mathematical Biology: Proceedings 



17 



of the Conference on Mathematical Biology and Dynamical Systems, World Scien- 
tific Publishing Co., 2007, pp. 15 - 38. 



[6] S. Schober, M. Bossert, Analysis of random boolean networks using the average 
sensitivity, ArXiv, \tt arXiv:nl.cg/0704.0197 (2007). 

[7] S. A. Kauffman, The Origins of Order: Self Organization and Selection in Evolution, 
Oxford University Press, 1993. 

[8] S. Kauffman, C. Peterson, B. Samuelsson, C. Troein, Random Boolean network 
models and the yeast transcriptional network, Proceedings of the National Academy 



of Sciences 100 (25) (2003) 14796-14799. ||doi : 10 . 1073/pnas . 20364291001 



[9] S. Kauffman, C. Peterson, B. Samuelsson, C. Troein, Genetic networks with can- 
alyzing Boolean rules are always stable, Proceedings of the National Academy 
of Sciences of the United States of America 101 (49) (2004) 17102-17107. 



doi: 10. 1073/pnas. 04077831011 



[10] T. P. Peixoto, The phase diagram of random boolean networks with nested canaliz- 
ing functions, The European Physical Journal B - Condensed Matter and Complex 
Systems 78 (2010) 187-192, 10.1140/epjb/e2010-10559-0. 

[11] C. H. Waddington, Canalization of Development and the Inheritance of Acquired 
Characters, Nature 150 (1942) 563-565. |doi : 10 . 1038/150563a0l 

[12] S. E. Harris, B. K. Sawhill, A. Wuensche, S. Kauffman, A model of transcriptional 
regulatory networks based on biases in the observed regulation rules, Complexity 



7 (4) (2002) 23-40. doi: 10. 1002/cplx. 10022, 



[13] F. Li, T. Long, Y. Lu, Q. Ouyang, C. Tang, The yeast cell-cycle network is robustly 
designed, Proceedings of the National Academy of Sciences of the United States of 
America 101 (14) (2004) 47814786, PMC387325. j doi : 10 . 1073/pnas . 0305937101| 

[14] M. I. Davidich, S. Bornholdt, Boolean network model predicts cell cycle sequence of 



fission yeast, PLoS ONE 3 (2) (2008) el672. |doi : 10 . 1371/ journal . pone . 0001672 



[15] J. G. Klotz, R. Feuer, K. Gottlieb, O. Sawodny, G. Sprenger, M. Bossert, M. Ederer, 
S. Schober, Properties of a Boolean network model of Escherichia coli, in: Proc. 
of the 8th International Workshop on Computational Systems Biology (WCSB), 
Zuerich, Switzerland, 2011. 



18 



[16] D. Murrugarra, R. Laubenbacher, Regulatory patterns in molecular in- 
teraction networks, Journal of theoretical biology 288 (2011) 66-72. 
doi:10.1016/j . jtbi . 2011 . 08 . 015 

[17] M. Gabbouj, P.-T. Yu, E. J. Coyle, Convergence behavior and root signal sets of 
stack filters, Circuits, Systems, and Signal Processing 11 (1) (1992) 171-193. 

[18] S. Cook, C. Dwork, R. Reischuk, Upper and lower time bounds for parallel random 
access machines without simultaneous writes, SIAM J. Comput. 15 (1) (1986) 87- 
97. |doi: 10 .1137/02150061 

[19] I. Benjamini, G. Kalai, O. Schramm, Noise sensitivity of boolean functions and 
applications to percolation, Publications mathematiques de 1'IHES 90 (1999) 5-43. 

[20] I. Shmulevich, Average sensitivity of typical monotone boolean functions (2005). 
larXiv : arXiv : math/0507030[ 

[21] Y. Li, J. O. Adeyeye, D. Murrugarra, B. Aguilar, R. Laubenbacher, Boolean nested 
canalizing functions: a comprehensive analysis (2012). arXiv: 1204.5203. 

[22] R. R. Bahadur, A representation of the joint distribution of responses to n dichoto- 
mous items, in: H. Solomon (Ed.), Studies on Item Analysis and Prediction, no. 6 
in Stanford mathematical studies in the social sciences, Stanford University Press, 
Stanford CA, 1961, pp. 158-176. 

[23] M. L. Furst, J. C. Jackson, S. W. Smith, Improved learning of AC0 functions, in: 
Proceedings of the Fourth Annual Workshop on Computational Learning Theory, 
Morgan Kaufmann Publishers Inc., Santa Cruz, California, United States, 1991, pp. 
317-325. 

[24] J. G. Klotz, S. Schober, M. Bossert, On the Predecessor Problem in Boolean Net- 
work Models of Regulatory Networks, in: Proc. of the 3rd International Conference 
on Bioinformatics and Computational Biology (BICoB), New Orleans, Louisiana, 
USA, 2011, pp. 67-73. 

[25] R. Heckel, S. Schober, M. Bossert, Harmonic analysis of boolean networks: Deter- 
minative power and perturbations, arXiv:1109.0807 (Sep. 2011). 

[26] M. Ben-Or, N. Linial, Collective coin flipping, robust voting schemes and minima 
of banzhaf values, in: FOCS, IEEE Computer Society, 1985, pp. 408-416. 



19 



[27] N. H. Bshouty, C. Tamon, On the fourier spectrum of monotone functions, Journal 
of the ACM 43 (4) (1996) 747-770. |doi : 10 . 1 145/234533 . 2 34564, 



20 



